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Abstract: In this paper we study the quantum wave packet of the non-linear Gross- 
Pitaeviskii equation. 

1. Introduction 

In the present work we investigate the quantum wave packet of the one-dimensional 
non-linear Gross-Pitaevskii equation with the potencial V(x, t) given by: 

V(x, t) = \ m uj 2 {t) x 2 , (1.1) 

which is the time dependent harmonic oscillator potential. 

2. Gross-Pitaeviskii Equation 

Em 1961 [1,2], E. P. Gross and, independently, L. P. Pitaevskii proposed a non-linear 
Schrodinger equation to represent time dependent physical systems, given by: 
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i h ^ii = _ J!L + I m uj\t) x 2 ^(x, t) + g\ ^{x, t) | 2 ^(x, f) , (2.1) 

where ip(x, t) is a wavefunction and g is a constant, constante. 

Writting the wavefuncition if>(x, t) in the polar form, defined by the Madelung-Bohm 
transformation [3, 4], we get: 

ip(x, t) = (f)(x, t) e i 5 ^ *) , (2.2) 

where S(x , t) is the classical action and <f>(x, t) will be defined in what follows. 

Substituting Eq.(2.2) into Eq.(2.1) and taking the real and imaginary parts of the 
resulting equation, we get [5]: 

t + = » (2-3) 

h f + \ m v U f ) + I ™ ^ + K« + ^gp = , (2.4) 
9j t + V + ^ = " ^ £ (K« + ^cp) , (2.5) 

where: 

p(x, £) = </> 2 (x, i) , (2.6) (quantum mass density) 
v qu (x, t) = A MgJl , (2.7) (quantum velocity) 

K«(x, i) = - ^ ^ = - 2~£^ , (2.8a,b) (Bohm quantum potential) 

and 

Vgp = 99- (2-9) (Gross-Pitaevskii potential) 
3. Quantum Wave Packet 

In 1909 [6], Einstein studied the black body radiation in thermodynamical equilib- 
rium with matter. Starting from Planck's equation, of 1900, of the radiation density and 
using the Fourier expansion technique to calculate its fluctuations, he showed that it exhibits, 
simultaneously, fluctuations which are characteristic of waves and particles. In 1916 [7], ana- 
lyzing again the black body Planckian radiation, Einstein proposed that an electromagnetic 
radiation with wavelenght A had a linear momentum p, given by the relation: 



p=$, (3.1) 
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where h is the Planck constant [8]. 

In works developed between 1923 and 1925 [9] de Broglie formulated his fundamental 
idea that the electron with mass m, in its atomic orbital motion with velocity v and linear 
momentum p = m v is guided by a "matter wave" (pilot-wave) with wavelenght is given 
by: 

A = \ ■ (3.2) 

In 1926 [10], Schrodinger proposed that the "pilot-wave de Brogliean" ought to obey 
a differential equation, today know as the famous Schrodinger's equation: 

i % | *(f, t) = H *(r, f) , (3.3a) 

where H is Hamiltonian operator definied by: 

H = ^ + V(r, t) , (p = -ikV), (3.3b,c) 

where V is the potential energy. In this same year of 1926 [11] Born interpreted the Schro- 
dinger wave function ^ as being an amplitude of probability. 

4. The Quantum Wave Packet of the non-linear Gross-Pitaeviskii equation 



Initially, let us calculate the quantum trajectory (x qu ) of the physical system repre- 
sented by the eq.(2.1). To do this, let us integrate the equation given by [12] [remember that 
q(t) = <*>]: 



where: 



and: 



v qu (x, t) = [x - q(t)\ + q(t) , (4.1) 



= h k(t)+ hoi . (4-2) 

cr(t) m yj TV a(t) 



k(t) = [ 2~~2j7\ H ] • (4.3) 



2 (t) 

So, remembering that J = in z, in (|) = in x — in y , and in x y = 
in x + in y], we have: 

Ms, t) = ^ = 4^ [x - q(t)} + q(t) , ^ 
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dx„ 



dt 



d[x qu (t) - q(t)] 

[Xqu(t) ~ <?(*)] 



&(t) dt da(t) 
2 cr(t) ~ 2 cr(t) 



ft d[x qu (t') - q(t')} 
J o [x qu {t>) - q(t')\ 



rt da(t') 
J O 2 cr(t') 



" J V [x qu (0) - g(0)] j - 2 Ul 



<r(0) 



<r(0) 



1/2 



g(«) 
<r(0) 



1/2 



[x gu (0) - q(0)} , (4.4) 



that represent the looked for quantum trajectory. 

To obtain the the quantum wave packet of the non-linear Gross-Pitaeviskii equation 
given by the eq.(2.2), let us expand the functions S(x, t), V(x, t) and V qu (x, t) around of 
q(t) up to second Taylor order [5]. In this way we have: 

S(x, t) = S[q(t), t] + S'[q(t), t] [x - q{t)\ + [x - g(t)] 2 , (4.5) 

V(x, t) = V[q(t), t] + V'[q(t), t] [x - q(t)} + ^MkA [x _ q[t) f > (4 . 6) 
V qu (x, t) = V qu [q(t), A + V^Mt), t] [x - q(t)} + ^SMhA [ x _ q{t )f , (4.7) 

V GP (x, t) = V GP [q(t), t] + V GP [q(t), t] [x - q(t)} + YM^iA [x _ q{t) f , (4 . 8 ) 

where (') and (") signifies, respectively, and J^. 

Differentiating the expression (4.5) in the variable x, multiplying the result by — , 
using the relations (2.6) and (4.1), and taking into account the polynomial identity property, 
we obtain: 



%_ dSjx, t) 
m dx 



m ( S'W% A + S"[q(t), t] [x - q(t)} ) 



= v qu (x, t) = [ [x - q(t)} + q(t) = - 

S'[q(t), A = ^ , S "[q(t), t] = f [ 4% ] . (4.9a,b) 
Substituting the expressions (4.9a,b) in the equation (4.5), results: 



S(x, t) = S (t) + 



m q(t) 



[x ~ q{t)] + fn 



a(t) j 



[x - q{t)f , (4.10) 



where: 



S (t) = S[q(t), t] , (4.11) 



is the quantum action. 



f = 0): 



Differentiating the eq.(4.10) in relation to the time t, we obtain (remembering that 



as 
at 



S (t) + | ( [x - q{t)] ) + | fk 



*(t) 



[x - q(t)f 



§ = So® + ^[X- q(t)} - HLfli + 

+ fk iW> - S® \? - ^ 2 *(*)] • ( 4 - 12 ) 

Considering that [12]: 



t) = ^p(x, t) = [7T (7(t)]- ^ e" 



2 <x(t) 



fx — x(t)] 

p(x, t) = [tt (r(i)]- 1/2 e ^) , (4.13) 



let us write V qu in terms of [x — q(t)]. Initially, we calculate the following differentiations: 
g = £ ( br .(()]" ■/« e - ^ ) = [, „(,)]- V« e - ^iL ( _ ) _ 



dx 



- [n a(t)}- V4 e - I^jM 



<T(t) 



l:L: JL ( _ [?r ff(t)] - 1/4 g - ^f- j^fll 



cte 2 cte 



<r(t) 



) = 



= - [tt a(t)}- V4 e 



2 o-(t) 



_9 / [s - 9(f)] \ 
9a; V o-(t) 7 



T(t) 



<r(t) 



Substituting the relation (4.14) in the equation (2.7b), taking into account the ex- 
pression (4.7), results: 



V qu {x, t) = V qu [q(t), A + V^Mt), A [x - q(t)} + V -^f^ [x - q{t)f 

*U*> *) = T^W) [* - m° - Tl&W) [* - (4-15) 
Comparing the eqs. (4.7) and (4.15), we have: 

t] = Y £ Wy V^ u [q(t), t] = 0; t] = -^Stfy (4.16a,b,c) 

Besides this the eq.(4.6) will be written, using the eq.(l.l) in the form: 
V(x, t) = V[q(t), t] + V'[q(t), *] [x - q(t)] + ™^ [x - q{t)f -> 

V(x, t) = \ m to 2 (t) q 2 (t) + 

+ ( m u\t) q(t) ) [x - q(t)} + f u\t) [x - q{t)f . (4.17) 

Comparing the eqs. (4.6) and (4.17), results: 

V[q(t), t] = ;m4) 9 2 (i); V'[q(t), t] = m u\t) q(t); (4.18a,b) 
V"[q(t), t] = mu 2 (t). (4.18c) 
Now, let us expand the eq. (4.13) around of q(t) to second Taylor order: 
p(x, t) = p[q(t), t] + p'[q(t), t] [x - q(t)\ + [x _ q{t)? _» 

p(x, t) = [n a(t)}- V2 e" + [* «r(t)]- ^ e * [x «f> ? + 

+ [* e" ( - ^ + ) [x - ,(,)]>. 

Considering only quadratics terms, results: 

p(x, t) = [tt t7(t)]" V* e" ( 1 + J^jQ£ ). (4.19) 

By using the eqs. (2.9), (4.8) and (4.19), we have: 
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V aP [q(t), t] = g[n a(t)]~ ^ e ( 1 + ) = 

= V GP [q(t), t] + V^ P [q{t\ t] [x - q(t)\ + YhMhA [x _ ,( t )]2 



K?p[g(*), *] = <? [tt a(t)]~ V2 e - V£ P [ ? (t), i] = 0; (4.20a,b) 



awl: 



V G p[g(t), t] = 2 3 ^ — 



(4.20c) 



Inserting the eqs. (2.7), (4.1; 4.5-8; 4.12; 4.14; 4.16a-c; 4.18a-c; 4.20a-c) into the 
eq.(2.4), we obtain, remembering that S (t), a(t) and q(t): 



h § + \ m v\ u + V + V qu + V GP 



= h 



So + V (* - «) " ^ + fn ( f " £ ) (* - ?) 5 



to tjr 
2 a 



( f ) (* - ?) 



+ | m 



2^ ) (z - 9) + 9 



+ 



+ \ m uj 2 q 2 + m u 2 q 2 (x — g) + y u 2 (x — q) 2 + 



2 m a 



2 m 



^ (x - q) 2 + g [tt ff (t)]- ^ e" l2L ^ L ( 1 + ) = . (4.21) 



Since (x — q)° = 1, we can gather together the above expression in potencies of 
(x — q), obtaining: 



h S — m q 2 + \ m q 2 + \ m u 2 q 2 + 2 — + <? [7r a] 1/2 e l ~ (x — g)° + 



1/2 P - 



+ 



m g — m q j-^ + m q ^ + m u z q 



[x 



- q) + 



m <t _ m (T 

4 cr 4 cr 2 



+ 



m <t 2 _i_ m uj 2 
8 a 2 ~T 2 



7i 2 , 2 g [tt 



1/2 ^t) 



2 m <r 2 



<T(t) 



(a; - g) 2 = . (4.22) 



As the above relation is an identically null polynomium, the coefficients of the po- 
tencies must be all equal to zero, that is: 



S (t) 



1 -2 1 2 2 

2 m g — ^ m u q 



2 m a 2 



- g[ 



1-1/2 - ^ - 9l " 

7T cr e 
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So® = \ 



1-2 1 2 2 Ti 

\mq 2 - \tuuj 1 q 2 - 



q + u 2 {t) q = , (4.24) 



9 P 



(4.23) 



2 cr 



a 2 
4 a 2 



+ UJ 2 + 



lJUL 



(4.25) 



Assuming that the following initial conditions are obeyed: 

q(0) = x , q(0) = v , cr(0) = a Q , <x(0) = b Q , (4.26a-d) 
and that [see eq.(2.7)]: 

S o (0) = ™J?^ ? (4.27) 
the integration of the expression (4.23) will be given by: 

S (t) = \ fl dt'\\m q 2 {t') - \ m u 2 {t') q 2 {t')- 



£m ~ 9 Pi*, t') 



2 m a(t') 



I m v x 
^ h 



(4.28) 



Taking into account the expressions (4.9a,b) and (4.28) in the equation (4.10) results: 



S(x, t) = \ J* dt'\\m q 2 (t>) - \ m u 2 (t>) q 2 (t>) 



2 m cr(t') 



9 p(x, t') 



+ 



+ 



m v x , m q(t) 



+ num[x- q(t)} + ?ft[z- q(t)} 2 . (4.29) 



This result obtained above permit us, finally to obtain the wave packet for of the 
non-linear Gross-Pitaeviskii equation. Indeed, considering the equations (2.2), (.6), (4.13) 
and (4.29), we get: 



V(x, t) = [?r a(t)\- 1/4 exp 



( m$ - ji® ) ^ - q w 



X 



x exp 



ij2 Y M [x - q{t)\ + 1 m 1° x ° 



X 



x exp 



if o dt'\\m q 2 (t>) - \ m u 2 (t') q 2 (t>) - 



2 m a(t') 



- 9 P{x, ?) 



(4.30) 



The eq. (4.30) we show that when q = 0, then we obtains the eq. (3.3.2.25) of the 
Reference [5], if a(t) = 2 a 2 (t), q(t) = X(t) and \ m u 2 (t) q 2 (t) = V[X(t)}. 
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